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ABSTRACT 

We derive the relative distances to the Virgo, Fornax, and Coma clusters of galaxies 
by applying the D n -a and the Fundamental Plane (FP) relations to the data of the 
homogeneous samples of early-type galaxies studied by Caon et al. (1990, 1994), Lucey 
J*"-. ' et al. (1991a, b), and J0rgensen et al. (1992, 1995a, b). The two distance indicators give 

■ consistent results, the relative distance moduli to Fornax and Coma with respect to 

Virgo being A^ipv = (0.45 ± 0.15) mag and A/icv — (3.55 ± 0.15) mag respectively. 
The formal error on Afi may be as small as 0.07 mag (~ 3% in distance), provided that 
all the sources of bias are taken into account and a correct statistical approach is used. 
Unfortunately, much of the actual uncertainty in the relative distance of the clusters 
(~ 12 — 15%), is due to the existence of systematic departures in the measurements of 

• the velocity dispersions among the various datasets, and to the corrections for aperture 

effects. 

The above result for the Fornax cluster is supported by the L—a—fj, relation and, 
£^ ■ with lesser accuracy, by the log(m)-log(r e ) relations. Our value of A/ipv is in fair 

^ \ agreement with the one derived using planetary nebulae and SNe-Ia, while is in open 

contrast with that coming from surface brightness fluctuations, globular clusters lumi- 
nosity function, and infrared Tully-Fishcr relation. In our data Coma appears slightly 
nearer than indicated by the other distance indicators, but now a better agreement 

■ with the Tully-Fisher relation seems to exist. 

5— i ' We show that for the galaxies of the Virgo and Fornax clusters the residuals of the 

D n -a relation do not correlate with the effective surface brightness (fi) c - There is also 
no correlation of the residuals of the D n -a and FP relations with the total luminosity 
of the galaxies, with the cllipticity and with the isophotal shape parameter 04. Instead, 
a correlation seems to exist with the maximum rotation velocity of the galaxies, V m , 
with the (V/a) ratio, and with the exponent m of the r x / m fit to the major axis light 
profiles of the galaxies. If confirmed, these effects introduce a systematic bias in both 
relations when used as distance indicators. 

Key words: Galaxies: clusters: Virgo, Fornax, and Coma - galaxies: distances - 
galaxies: elliptical and lenticular - galaxies: kinematics and dynamics 



1 INTRODUCTION 

Empirical scaling laws built using the main structural pa- 
rameters of galaxies such as characteristic radii, mean sur- 
face brightnesses, and central velocity dispersions, have be- 
come quite fashionable in the last decade. In particular, the 
relations holding for early-type galaxies are believed to be 
useful tools for the understanding of the structure, forma- 



tion, and evolution of the kinematically hot stellar systems 
(hereafter early-type galaxies) , and to be also quite accurate 
distance indicators (DIs) for these objects. 

One of such relations is the so-called Fundamental 
Plane (hereafter FP; Djorgovski & Davis 1987), which, in 
its original formulation, binds the central velocity disper- 
sion, <7o, and the average surface brightness within the ef- 
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fective isophote (that encircling half of the total light of 
a galaxy in a given photometric band), (/Lt) e , to the mean 
radius of this same isophote, named effective radius r c : 
log(r e ) = a log(ao) + b{fj,) e + c. 

A second scaling law is the D n -a relation (Dressier et 
al. 1987): D n oc (Tq, where D n is the diameter of the circular 
aperture within which the average corrected surface bright- 
ness equals a fixed value (fj,) n usually 20.75 B mag arcsec - , 
and 7 is a constant of the order of 3/2. To a first order ap- 
proximation the D n -a relation is an edge-on view of the 
FP (Dressier et al. 1987) . In this case the surface brightness 
parameter is hidden in the definition of D n . 

The slope and the small thickness measured for the FP 
and the D n -a relations imply a narrow range of density 
structure, velocity anisotropy, and mean age among ellipti- 
cal galaxies, which are excellent properties for DIs. Actually 
these relations have provided much of the present informa- 
tion on the large scale flows in the local universe (Lynden- 
Bell et al. 1988). 

Two complementary strategies have been developed. 
On the one hand the enormous progress in the observa- 
tional techniques and data reduction procedures has pro- 
duced larger and larger samples and increasingly better dis- 
tance determinations of galaxies. On the other hand much 
effort has been put in the last years in quantifying internal 
errors, biases, calibration uncertainties, second parameter 
corrections, environmental influences, and behaviour of the 
residuals in both relations; see e.g. Lucey et al. (1991a; b 
= LGCT), Jacoby et al. (1992), j0rgensen et al. (1992; 
JFK 1995a,b), Davies et al. (1993), Saglia et al. (1993), van 
Albada et al. (1993). In conclusion, from the original 32% 
uncertainty given by the Faber & Jackson (1976; FJ) rela- 
tion, the error has been progressively reduced by the coming 
of the D n -a and of the FP relations, being nowadays as low 
as 10% in the relative distances between clusters (Jacoby et 
al. 1992). 

Despite this success, a number of problems are still 
open. Which of the two relations, D n -a and FP, is more 
suiteable as DI ? Or, in other words, which of the two is less 
affected by observational errors and biases ? Can SO galaxies 
be included in these relations ? Many papers have addressed 
the above questions with contradictory results. 

JFK state that the FP should be preferred, because the 
D n -a relation has a larger scatter and suffers from a surface 
brightness bias. Van Albada et al. (1993) find that the D n -a 
relation has the same accuracy of the FP, provided that a 
proper surface brightness correction is applied. Lucey et al. 
(1991a) and LGCT derive contradictory evidences for such a 
bias in their samples, and Saglia et al. (1993) show that the 
scatter around both relations can be reduced by rejecting E 
galaxies with inner disk components. Finally, Davies et al. 
(1993) point out that only galaxies in a restricted range of 
surface brightness should be used in the D n -a relation in 
order to avoid contaminations. 

We argue here that the D n ~a and the FP relations have 
substantially the same accuracy - at least as long as the 
distance to nearby clusters is concerned - by applying them 
to the measure of the relative distances to the Fornax and 
Coma clusters with respect to Virgo. The motivation of this 
work is twofold. Firstly, the relative distances to these clus- 
ters are far from being accurately known (cf . Section ^|) , and 
secondly, today we can exploit an accurate and homogeneous 



data-set of photometric and spectroscopic measurements for 
the galaxies of these three clusters: Caon et al. (1990, 1994, 
hereafter C 2 D), D'Onofrio et al. (1995), LGCT, JFK, and 
j0rgensen et al. (1995a, b). The study of a volume-limited 
sample of galaxies has the advantage that one can look at 
the intrinsic properties of the D n -a and FP relations by 
means of objects which are probably coeval, and reduce as 
much as possible the principal sources of errors, in particular 
those related to the unknown distances of the galaxies. 

The paper is organized as follows. Section ^| reviews the 
current determinations of the relative distances to the Virgo, 
Fornax, and Coma clusters, as measured by the relevant DIs. 
Our determinations of the relative distances to the three 
clusters using the D n -a and FP relations are presented in 
Sections ^| and ^ respectively, together with a discussion of 
the errors intrinsic to both methods and a comparison of the 
results with other DIs. In Section ^ we analyze the behaviour 
of the residuals of the D n -a and FP relations by exploiting 
the advantage of our volume-limited, accurate, and homo- 
geneous (in the methodological approach and in the data 
analysis) sample of galaxies, and we perform some further 
tests on the D n -a relation by varying the definitions of D n 
and a. In Section ^| we present the log(m)-log(r e ) relation 
used as DI for the Fornax cluster. Conclusions are drawn in 
Section 71 



2 THE RELATIVE DISTANCES TO THE 

VIRGO, FORNAX, AND COMA CLUSTERS 

In his review of the extragalactic distance scale, de Vaucou- 
leurs (1993) showed that the differential distance modulus 
(A[i) between the Fornax and the Virgo cluster is, on aver- 
age, A/xfv = +0.16 and +0.38 mag for the long and for the 
short distance scale respectively. In other words, according 
to this study, both scales place Virgo closer than Fornax. 

The trend is not confirmed by all DIs. Table 1 lists the 
measurements of A/ipv— /if — [iv, sorted according to the 
method used; we have essentially ignored the old measure- 
ments, already revisited by de Vaucouleurs. A quick look to 
the Table shows that the inner ring diameters of spirals, 
the globular cluster luminosity function (GCLF), the sur- 
face brightness fluctuations (SBF), the L-a-fj, relation, the 
IR Tully-Fisher relation (IR-TF), and the scale length of 
dE galaxies, place the Fornax cluster closer than Virgo, in 
marked contrast with the results from the brightest cluster 
members (BCM), the D n -a relation, the planetary nebulae 
luminosity function (PNLF), and the SNe-Ia. Taking a plain 
average, the A/ipy turns into a salomonical ~ 0.0, with a 
r.m.s. of 0.24 mag. 

A few comments to Table 1 are in order. The SBF 
method gives distances systematically lower by 10-30% than 
the D n -a relation, probably due to the adopted calibration 
(Sandage & Tammann 1995). The values of A/x obtained 
from the analysis of the GCLF strongly depend on the mor- 
phological types of the selected galaxies and on metallic- 
ity effects, while those provided by the PNLF method are 
in good agreement, within the errors, with the D n -a rela- 
tion (McMillan et al. 1993). For both methods however, the 
number of galaxies in common with our D n -a relation is 
quite small, and a conclusive comparison cannot be made. 
The filamentary and complex structure of the Virgo cluster 
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Table 1 : Distance modulus of the Fornax cluster relative to that of Virgo 





A i. 




Mctnoa 


Reference 


IN OtCS 


+0 


.48 ±0 


.20 


BCM lum. 


de Vaucouleurs 1977 


Uncertain estimate of the error 


-0 


.37 ±0. 


.30 


Ring diam. 


Buta & de Vaucouleurs 1983 


Based on inner rings of 4 spirals in Fornax and 17 in Virgo 


+0. 


.15 ±0. 


28 


D n -a 


Davies et al. 1993 


Calibrated through the Leo Group set at a distance of 10 Mpc 


+0. 


,13 ±0. 


20 


D n -o 


Dressier ct al. 1987 


The estimated error concerns galaxies in clusters 


+0. 


.25 ±0. 


.31 


D n -a 


Fabcr et al. 1989 


See McMillan ct al. (1993) 


-0. 


.14 ±0. 


.28 


L-a-fi 


Pierce 1989 


Assuming Leo Group at 10 Mpc 


+0. 


00 ±0. 


22 


L-T, 


Ferguson & Sandagc 1988 


See McMillan ct al. (1993) 


-0. 


,15 ±0. 


.13 


SBF 


Tonry 1991 


Assuming fj, = 24.43 mag for M31 and M32 and no correction 












for absorption 


-0. 


.17 ±0. 


.34 


SBF 


Buzzoni 1993 


Assuming M32 at 0.74 Mpc 


+0. 


.08 ±0. 


.11 


SBF 


Tonry 1991 


See McMillan ct al. (1993) 


-0. 


50 ±0. 


20 


GCLF 


Geisler & Forte 1990 


Match of the LF of NGC 1399 and M87 


-0. 


.11 ±0. 


,14 


GCLF 


Bridges et al. 1991 


Using the LF of NGC 1399 and a composite Virgo LF 


-0. 


.16 ±0. 


.15 


dE s. length 


Bothun ct al. 1989 


From the scale length of ~ 25 dwarf ellipticals in both clusters 


-0. 


25 ±0. 


23 


IR-TF 


Aaronson et al. 1989 


H-band TF for 16 spirals in Virgo and 7 in Fornax 


+0. 


.24 ±0. 


.10 


PNLF 


McMillan ct al. 1993 


From the PNLF of NGC 1316, NGC 1399, and NGC 1404 


+0. 


03 ± 0. 


10 


PNLF 


Ciardullo et al. 1991 


iFrom de Vaucouleurs (1993) 


+0. 


.15 ±0. 


.94 


SNe-Ia 


Hamuy et al. 1991 


B-band maximum magnitude of SNc 1980N and 1981D in 












NGC 1316 and SN 1984A in NGC 4419 


+0. 


26 ±0. 


71 


SNe-Ia 


Hamuy ct al. 1991 


Maximum V-band magnitude of the same objects 


+0. 


.09 ±0. 


.14 


SNe-Ia 


Hamuy et al. 1991 


Maximum H-band magnitude of the same objects 



probably accounts for the smaller distance found by the IR- 
TF relation, since the spiral galaxies have a different spatial 
distribution than ellipticals. The negative Afi value given 
by Pierce (1989) through the L-a-fi relation is probably 
the consequence of the small size of his sample: we obtain 
A/kfv = by applying the D n -a relation to the galaxies of 
our sample in common with his. This is a clear demonstra- 
tion that a "representative" statistical sample of galaxies in 
a cluster, selected according to well defined criteria, must be 
used to gauge relative distances through the D n -a and FP 
relations (see next section). 

In a similar way Table 2 reviews the determinations 
of the relative distance between the Virgo and the Coma 
clusters. The average differential distance mudulus is now 
3.75 ±0.18 mag, but the individual determinations of A/icv 
span over an interval of ~ 0.6 mag. Most of the measure- 
ments range however between 3.60 and 3.80 mag. If the cen- 
tre of the Virgo cluster is taken at 18.3 Mpc (Capaccioli et 
al. 1990), considering the broad range of values in Table 2, 
the distance to Coma ranges from 88.4 and 119.3 Mpc. 

We note that the D n -a relation of Dressier et al. (1987) 
gives results in good agreement with the SNe-Ia method and 
with the TF determinations. It is also apparent from the 
above Tables that the principal DIs provide quite different 
distances for each cluster. We will not try to give here an 
explanation for such discrepancies, since this would require 
a large number of distance determinations in common be- 
tween the various DIs. In the following we will address only 
the problems encountered with the D n -a and FP relations, 
when they are used as DIs. 



3 THE RELATIVE DISTANCES TO THE 

VIRGO, FORNAX, AND COMA CLUSTERS 
THROUGH THE D n -a AND THE FP 
RELATIONS 

There are various sources of errors affecting the determina- 
tions of the relative distances to nearby galaxy clusters with 
the D n -a and FP relations. The selection criteria of the 
galaxy sample play an important part in this game, together 
with the accuracy and homogeneity of the photometric and 
spectroscopic data, the definition of the variables involved 
in the relations, and the statistical regression analysis that 
is used. 

A number of questions still awaits for a clarification. 
Does it exist a "representative" sample of early-type galax- 
ies in clusters providing unbiased distance determinations ? 
Can SO galaxies be included in the D n -a and FP rela- 
tions ? Which are the "best" parameters to use ? For in- 
stance, shall we use the effective radius r e derived from 
a model-dependent fit to the light profiles (e.g. r 1 / 4 law), 
rather than the radius of the isophote which encircles half 
the total luminosity, whatever the light profile be? Shall 
we choose the central velocity dispersion <j or an average 
within a given radial range ? Is it better suited the diameter 
D n at fiB = 20.75 mag arcsec _2 or the diameter relative to 
a different mean surface brightness level ? We have no clear 
answers to any of such questions yet. 

Unfortunately, the accuracy and homogeneity of the 
spectroscopic data is still low compared with those of pho- 
tometric data. Extended rotation curves and velocity dis- 
persion profiles do exist for a few galaxies only (usually the 
brightest ones), and in most cases the full velocity field of 
the objects is poorly known. Systematic and/or random dif- 
ferences in the velocity dispersion measurements are eas- 
ily found in the literature for the same objects and likely 
arise from the use of different methods of data reduction 
(Fourier Quotient, Cross Correlation, Fourier Correlation 
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Table 2: Distance modulus of the Coma cluster relative to that of Virgo 



A^cv 


Method 


Reference 


Notes 


+6.72 ± 0.30 


TV T £ r C\T 

Max. of 5 bNe 


dc Vaucoulcurs 1985 


Irom lablc 16B in dc Vaucoulcurs (1993) 


+3.65 ± 0.20 


D n -cr 


Dressier et al. 1987 


.brom a sample of 20 bs in Virgo and 28 fts m Coma 


+3.99 ± o.do 


D n -Mg 2 


~T"\ ........ 1 _ , . . 1 1 { \ i) "~7 

Dressier et al. 1987 


The error is the 35% r.m.s. uncertainty of the method 


+3.82 ± 0.10 


D„-IT 


leigelson & Babu 1992 


Data Irom LGC1 and Dressier ct al. (1987) 


+3.75 ± 0.20 


SNe-Ia 


Capaccioli et al. 1990 


Data Irom the Asiago biN Catalogue ol Barbon ct al. (1989) 


+3.60 ± 0.30 


SNe-Ia 


Capaccioli et al. 1990 


Using SNe— la in 4 E galaxies 


+3.69 ± 0.12 


IR-TF 


A f> TV T 111 / \ / ■ 

Aaronson & Mould 1986 


bee Tabic 4 of Capaccioli ct al. (1990) 


1 1 f?f~* 1 C\ "1 A 

+3.66 ± 0.14 


U — B of E 


bandage 1972 


Colour-magnitude (CM) diagram. Revised to 3.5+0.2 by 








Aaronson ct al. (lyol) 


+3.70 ± 0.17 


B and 11 11 


Giraud 1986 


TT 1 " 1 J? J? J_l TITT* 1 j_ " 

Hybrid lorm ol the 1 1 relation 


+3.76 ± 0.12 


L-cr 


Lucey 1986 


through hxed metric apertures 


+3.75 ± 0.18 


V26-C 


Dressier 1984 


V26 Irom bandage & Visvanathan (1978) 


+4.U / ± 0. / 


Reduced radii 


ouacnus tyyi 


Mean error for a bright isolated galaxy 


+3.42 ± 0.22 


B-band TF 


Rood & Williams 1992 


For M31— like galaxies 


+3.70 + 0.09 


Mean of 3 


Bower ct al. 1992 


Mean of CM diagram, eolour-diamcter and FJ relations 


+3.59 + 0.06 


CM diagr. for E 


Bower ct al. 1992 


As above for E galaxies only 


+3.80 + 0.06 


Mult, analysis 


Vader 1986 


Comparison between the mass-luminosity relations 


+3.60 + 0.17 


IR CM diagram 


Christcnscn 1991 


Unbiased value for spiral galaxies 


+3.90 + 0.10 


Mean of 9 


Tammann & Sandagc 1985 


Mean of 9 determinations from the literature 


+4.18 + 0.19 


Mean of 3 


Sandage & Tammann 1982 


bee Table 4 of Capaccioli ct al. (1990) 


+3.60 + 0.40 


B-band TF 


Fukugita ct al. 1993 


Virgo centre placed at 15 Mpc 



Quotient, etc.) and to the corrections for the aperture size 
effects (j0rgensen et al. 1995b). 

Presently the Virgo, Fornax, and Coma clusters offer 
the best studied galaxy samples for addressing the above 
questions, since the available data approximately satisfy the 
required criteria of homogeneity and accuracy needed for the 
analysis of these problems. 

3.1 The samples 

The C 2 D set of B-band photometric parameters for Virgo 
E and non-barred SO galaxies is 80% complete down to 
Bt = 14 mag (according to the catalog of Binggeli et al., 
1985, hereafter BST), and that of Fornax is 100% complete 
down to Bt = 15 mag (according to the membership list of 
Ferguson, 1989). Assuming the same distance of 18.3 Mpc 
to both clusters, the samples are respectively complete (and 
volume- limited) down to Mb = —17.3 and Mb = —16.3 
mag. 

The D n diameters, effective radii r c , and effective mean 
surface brightnesses (/i) e , are given by C 2 D with an accuracy 
of about 5%, 20%, and ~ 0.3 mag (~ 30%) respectively. The 
above parameters have not been corrected for seeing, galac- 
tic absorption, internal extinction and redshift. All these 
corrections are very small for the galaxies in both clusters, 
as explained in C 2 D. 

For the photometric parameters of the Coma galaxies 
we made use of the data of LGCT, JFK, and J0rgensen et 
al. (1995a), which have a similar degree of homogeneity and 
accuracy as C 2 D. These data have been corrected by the 
authors for the above cited effects. 

Unfortunately, the spectroscopic data are far less homo- 
geneous. We have measured the rotation curves V(r) and 
velocity dispersion profiles <r(r) of 15 galaxies of the C 2 D 
Fornax sample (D'Onofrio et al. 1995), while for the Virgo 
galaxies we collected central velocity dispersion measure- 
ments from the literature: Whitmore et al. (1985, WMT), 
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Figure 1. Panel a) Distribution of the kinematical subset of the 
Virgo galaxies of the C 2 D sample in the logD„-logcr plane. The 
solid line gives the fit to the "best sample" , obtained by rejecting 
the objetcs deviating largely from the mean relation. Panel b ) The 
same diagram for the Virgo ("best sample") and Fornax early- 
type galaxies. The dashed line is the fit of minimum scatter for 
the galaxies of the Fornax cluster, adopting the same slope of the 
Virgo relation. 



Dalle Ore et al. (1991, DFJS), Davoust et al. (1985, DPV), 
Faber et al. (1989, 7-Sam), McElroy (1995). For the Virgo 
cluster we were able to find kinematical data for 12 galax- 
ies only (10 E and 2 SO): NGC4374, 4486, 4636 (Davies 
& Birkinshaw 1988), NGC4365, 4621 (Bender 1988a,b), 
NGC4387, 4478 (Davies et al. 1983), NGC4406, 4472, 
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4649 (Franx et al. 1989), NGC4459 (Peterson 1978) and 
NGC4473 (Young et al. 1978). 

For the galaxies of the Coma cluster we used the cen- 
tral velocity dispersions quoted by LGCT, JFK (who homo- 
geneized the data of Faber et al. 1989 and Dressier 1987), 
and J0rgensen et al. (1995b, who increased the JFK sample 
and obtained average velocity dispersions from a fixed and 
a variable, within r c /4, aperture of the slit). 

The accuracy of the kinematical data is more difficult to 
evaluate. Much of the offset in the distance determinations 
are likely caused by systematic differences in the a mea- 
surements. Such differences are difficult to evaluate due to 
the heterogenous nature of the samples and of the methods 
adopted to gauge the central velocity dispersion. 

According to Jacoby et al. (1992) the errors on a likely 
range between 6 and 14%. ft is clear that such an estimate 
does not account for possible local deviations from the over- 
all dynamical behaviour due for instance to singularities, 
which may cause the appearance of central spikes in the ve- 
locity dispersion profile (which behave as systematic errors 
on a). 

Of the 52 early-type galaxies in the Virgo sample (29 
SO and 23 E), 13 SOs and 4 Es have no cr , 3 SOs have 
(To < 100 km s , and 4 SOs and 2 Es have anomalous 
high/low redshifts with respect to the mean of the cluster 
(czo = 1179 ± 17, Sandage & Tammann 1995). One galaxy 
(NGC 4370) has been eliminated because its central surface 
brightness is fainter than /is = 20.75 (implying D n = 0). 
In Fornax we have 12 SOs and 16 Es: 7 SOs and 7 Es have 
no measurements of oo, and 2 SOs have uq < 100 km s . 
There are no peculiar redshifts. 

3.2 The D n —a relation for the Virgo and Fornax 
clusters 

Fig. |l| plots the D n -a relation for the C 2 D sample. The 
upper panel shows only the Virgo galaxies, with different 
symbols according to their morphology: open circles for Es 
and triangles for SOs. Here the central velocity dispersions of 
the Virgo galaxies are the plain averages of the data found 
in the literature. 

It is apparent that E and SO galaxies follow the same 
relation. There are, however, 8 Virgo galaxies that largely 
deviate from the mean relation (indicated by the solid line; 
see below). 

The discrepant galaxies have peculiar redshifts and/or 
high inclination to the line of sight. Three of them have <to < 
100 km s -1 and, following the recipe of the 7-Sam have been 
discarded. Of the 9 SOs that follow the mean relation, 3 have 
apparent flattening higher than S0(5) (according to the BST 
classification) and one of them has a peculiar redshift. Two 
galaxies fitting the mean relation are rounder than S0(5) but 
have discrepant redshifts. Therefore it seems that inclination 
and redshift do not determine univocally the galaxies that 
follow the D n -a relation and are not useful selection criteria 
for establishing a "representative" sample. 

Table 3a presents the morphological type, the axis ratio, 
the recession velocity and the BST membership code for the 
Virgo galaxies which deviate significantly from the mean 
D n -a relation, and Table 3b the same data for those SOs 
that instead fit the relation. In the following we will reject 
the 8 discrepant galaxies when computing the D n -a and FP 



Table 3a: Virgo galaxies that do not fit the D n -cr relation 



Id cut 


1V11_/1 ]_J11. 

Type 


(-) 

\aJ 25 


[ km s 1 ] 


VI r»rn li/irth i t i 

1V1L1I1LJ\1.1 Gill YJ 


NOC 41 fi8 

li VJ*J rrlUO 


E2 


0.79 


2316 




NGC 4200 


S0(4) 


0.61 


2376 




NGC 4261 


E2 


0.80 


2200 


- 


NGC 4270 


S0(6) 


0.46 


2347 




NGC 4281 


S0(6) 


0.56 


2711 




NGC 4342 


S0(7) 


0.54 


714 




NGC 4417 


S0(7) 


0.62 


832 


M 


NGC 4550 


S0(7) 


0.31 


378 


M 


Table 3b: Virgo 


SOs fitting the D n —(T 


relation 


Went. 


Morph. 


(") 


cz 


1V1^111U\_.1 Olll l-J 




Type 


\aj 25 


[ km s _1 ] 




NGC 4339 


S0(0) 


0.86 


1287 




NGC 4377 


S0(3) 


0.85 


1371 


M 


NGC 4459 


S0(2) 


0.84 


1210 


M 


NGC 4526 


S0(6) 


0.37 


533 


M 


NGC 4552 


S0(0) 


0.85 


321 


M 


NGC 4570 


S0(7) 


0.31 


1730 


M 


NGC 4578 


S0(4) 


0.70 


2284 


M 


NGC 4638 


S0(7) 


0.70 


1147 


M 


NGC 4649 


S0(2) 


0.81 


1144 


M 



Notes to Tables 3a, b: Morphological types, recession velocities 
and membership are from Binggeli et al. (1985); axis ratios at 
Hb = 25 are from C 2 D. 



relations (the same galaxies, with the possible exception of 
NGC 4200, also show a large scatter with respect to the FP); 
they are likely foreground and background objects. 

The solid line in Fig. |l| has been computed through 
a linear least-square fit of the best Virgo sample, ignoring 
errors on both coordinates. We get: 

\og(D n ) = (1.26 ± 0.06) log(o-o) - (1.11 ± 0.14) (1) 

with a r.m.s. of 0.06 in log(D n ), corresponding to a 15% 
uncertainty in the average distance. 

With a more robust fit which takes into account the 
errors on both coordinates (Fasano & Vio 1988), it is: 

log(X)„) = (1.31 ± 0.07) log(ao) - (1.24 ± 0.16) (2) 

and the same r.m.s. fit as aboveQ. 

Taking into account the above quoted errors, the intrin- 
sic scatter of the D n -a relation is 0.05, corresponding to a 
distance uncertainty of ~ 12% per galaxy. 

In the lower panel of Fig. [l| the Virgo and Fornax galax- 
ies are indicated with different symbols. The solid line is our 
eq. §. 

If we assume that the D n -a relation for the Fornax 



* Note that the errors on both coordinates (Alog(Dn) ~ 0.02 
and Alog(cr) ~ 0.04) are approximately equal, and the use of a 
least square fit that weights both variables does not produce a 
significantly different result. 
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cluster has the same slope as that for Virgo (a reason- 
able approximation if the environment does not influence 
the D„-a relation), the fit providing the minimum scat- 
ter for the Fornax galaxies gives Alog(D n ) = 0.09, and 
a ratio -D„(Virgo)/D n (Fornax) = 1.23, corresponding to 
A/^fv = 0.45 mag. 

This result disagrees with Davies et al. (1993), who 
found D n (Virgo)/D n (Fornax) = 1.07, and a A/iFV = 0.15 
mag. The origin of such a discrepancy (15% in distance) 
is not in the values of D„, since the two samples compare 
quite well, but rather in the values of a. For 7 of the 8 Fornax 
galaxies we have in common, the values of ao are larger in 
our sample. The discrepancy between the spectroscopic data 
of the two samples increases for the objects with ao > 150 
km s" 1 (D'Onofrio et al. 1995). 

The non-homogeneity of the velocity dispersion data is 
also clear from Fig. [I]. Note that the scatter for the Virgo 
objects is lower than for Fornax: the r.m.s. of the least square 
fit of the two clusters is 0.06 for Virgo and 0.13 for FornaxQ 

This is at variance with the expectation based on the 
morphology of the two clusters: diffuse (Virgo) and com- 
pact (Fornax), according to BST and Ferguson (1989). If 
one assumes that the two clusters have spherical symmetry, 
an intrinsic ~ 7% and ~ 11% scatter is expected for Fornax 
and Virgo respectively, according to their angular sizes. 

One might speculate that the small scatter of the Virgo 
data is due to the averaging of several measurements of a, 
obtained by different authors using different methods of data 
analysis on spectra taken with different slit apertures. What 
does it happen if we do not average the central velocity dis- 
persion measurements ? This is shown in Fig. ^| where, for 
each source of ao data, we report the slope, intercept, and 
r.m.s. scatter of the best fit for the corresponding Virgo 
galaxies and the resulting A^ifv computed adopting the 
slope derived for Virgo. 

This experiment tells us two things. The first is that 
the small scatter of the Virgo data is independent of the 
databases (each catalog of velocity dispersions is made of 
averaged and homogenized data). If we use, however, the 
12 galaxies of the sub-sample with the best observed rota- 
tion curves and velocity dispersion profiles, the r.m.s. scatter 
from the best fit increases to 0.08 (20%). The second is that 
the relative distance moduli between the two clusters span 
a large interval, varying with the properties of the adopted 
Virgo sample. The null value of A/ifv found in our compar- 
ison with Pierce (1989) demonstrates this effect. 

Since most of the Virgo data come from the 7-Sam sam- 
ple, we decided to concentrate on the errors introduced by 
the aperture sizes on the relative distance between Virgo 
and Fornax, adopting the central velocity dispersions of their 
Lick and LCO-HI measurements (4" x l".5 and 4" x 4" slit 
aperture respectively). The two data-sets for precisely the 
same galaxies, give A/ifv = 0.35 and 0.55 respectively. This 
means that an accurate distance between the clusters could 
be obtained only after an homogenous comparison of the 
a measurements, corrected for the aperture size effects. A 
more quantitative idea of such corrections has been provided 
by J0rgensen et al. (1995b). 

t Excluding NGC 1316, the galaxy that has the largest discrepant 
value of (To, the r.m.s. lowers to 0.10. 



Finally, we calculate the distance between Fornax and 
Virgo through the L — a relation using the C 2 D B- 
magnitudes for the galaxies of the two clusters and the ve- 
locity dispersions from the catalog of McElroy (1995). We 
obtained A/ifv = 0.43 mag, in good agreement with the 
D n -a relation. 

3.3 The D n —a relation for the Coma cluster 

We used the data given by LGCT, JFK and J0rgensen et 
al. (1995a, b) to derive the distance to the Coma cluster. 
Unfortunately, the D n diameters measured by LGCT are in 
the V band (instead of B), at the mean surface brightness 
level of 19.8 mag arcsec -2 . However, since the average B — V 
colour of elliptical galaxies is approximately constant and 
equal to 0.95 mag, we can equally compare the two data 
sets. The same happens for the D n values of J0rgensen et 
al. (1995a) which are in the r band. In that case, using the 
photometric parameters in the B and r bands listed in JFK, 
we derived the relations: log(r^) = 0.968 log«) + 0.017 
and log(r^) = 1.006 log(r£) +0.026. The new parameters of 
j0rgensen, together with those given by LGCT, have been 
corrected by the authors for cosmological and seeing effects. 

For the LGCT and JFK samples the line of minimum 
scatter with respect to eq. |l] (Fig.g) gives Alog(D n ) = 0.68 
and Alog(D„) = 0.69 respectively, corresponding to a 
A/xcv of 3.4 and 3.45 mag. Considering a 5% aperture cor- 
rection to the observed ao passing from the Virgo and Coma 
galaxies (as in Davies et al. 1987), A/icv turns out to be 
~ 3.55 mag. We note that the use of the LGCT seeing cor- 
rected parameters does not influence very much the final 
result. 

We also calculated the relative distance between Coma 
and Virgo using the photometric parameters and the Coma 
velocity dispersions quoted by j0rgensen et al. (1995b), who 
give the average a within a fixed aperture of 3". 4 and within 
a variable aperture r e /4, which is connected to the physical 
size of the galaxies. It was impossible for us to get the same 
variables for the Virgo galaxies, since this requires the ve- 
locity dispersion profiles. We found A^tcv = 3.45 and 3.55 
mag respectively in the two cases, confirming the 12-15% 
distance uncertainty due to aperture effects. Since the r e /4 
aperture is generally lower than the 3". 4 slit for the Coma 
galaxies, we conclude that, by enlarging the aperture, one 
obtains a lower average value of a and consequently a lower 
differential distance modulus. 

It is also interesting to compare the D n -a results with 
the V26 — o relation of Dressier (1984, hereafter AD), who 
found A/icv = 3.75 mag. His result is particularly important 
here because Dressier considered two different slit apertures 
for the galaxies in Virgo and Coma in order to reduce as 
much as possible the aperture effect. The ~ 0.25 mag dif- 
ference of the two methods can be explained either by a 
systematic difference in D n (15%) or in ao (10%). 

We first compared the magnitudes of the Virgo galaxies 
in common with Dressier (he used the data of Sandage & 
Visvanathan 1978). ^,From the growth curves of the C 2 D 
sample, we derived the B26 magnitudes of our galaxies, and 
compared them with Dressler's data by transforming his V26 
magnitudes in the B— band. We adopted the color equation 
of Sandage & Visvanathan (1978): B — V = 1.25(6 — V) + 
0.22. We obtained B our = 1.03 B Dre - 0.80. Being only a 
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Figure 2. The D n -o relation using individual sources of <tq for the Virgo galaxies. Open and filled symbols indicate Virgo and Fornax 
objects respectively. The solid lines are the fit the Virgo data, for which we report intercept a, slope b and r.m.s.. The dotted line is our 
eq. hi The dashed lines are the fit for Fornax, assuming the same slope of Virgo, from which Afj, is computed. 



zero point offset in the B magnitudes, the resulting B26 — er 
relation for the Virgo- Coma galaxies produces a Apcv of 
3.80. Unfortunately, Dressier does not report the D n of his 
galaxies, and so a definitive comparison cannot be made. 
However, when we compared our data with Dressler's et al. 
(1987), we found a good agreement. This is an indication 
that the observed discrepancy is not due to a systematic 
variation of the photometric parameters. 

The discrepancy between the D n -a and the L — a re- 
lations can also be due to a systematic ~ 10% offset in the 
velocity dispersions. Part of this offset (Alog(c) = 0.028) 
has been found by LGCT in their comparison with the data 
of Dressier (unfortunately however, such a correction to the 
a values of Dressier brings the result of the B — a relation 
to A/xcv = 4). By comparing the central velocity disper- 
sions of the Coma galaxies from LGCT, Dressier (1984), and 
McElroy (1995, hereafter McE), we found A log cro(McE - 

LGCT) 0.017 and Alogao(McE-AD) ~ 0.014. For the 

Virgo galaxies we made the comparisons between Dressier 
(1984), McElroy (1995) and Davoust et al. (1985, here- 
fater DPV) obtaining Alog a (McE - AD) ~ 0.015 and 
Aloga (DPV - AD) ~ 0.04 (note that the DPV data have 
been homogeneized by the authors taking into account aper- 



ture effects, method of reduction, etc.). This means that the 
aperture effects can account for the observed difference. The 
large differential distance modulus is likely the result of the 
low central velocity dispersion that Dressier got for the Virgo 
galaxies. 

Looking at the data of the 7-Sam and of DPV, we dis- 
covered that in the 70 ™ 75% of the cases the use of a small 
aperture (slit width over face on corrected diameter D(0) 
less than 0.02) produces an higher central velocity disper- 
sion (01/ a s = 0.85), where the subscripts here stay for "large 
aperture" and "small aperture". This effect can be easily 
responsible of the scatter observed in A/icv, since it gives 
~15% difference in the central velocity dispersions. 

A further possible explanation of the difference between 
the L — a and D n -a relation resides in the surface brightness 
bias which affects the L — a relation. To test this effect we 
corrected the total B magnitudes of the galaxies adopting 
the relation AB = — 0.3{[i) e + 6.49, which represents the 
trend of the residuals as a function of the effective surface 
brightness in our data. The result is that we obtain a change 
in the slope and in the scatter of the B — a relation, but the 
differential modulus is not affected. 

The distance to Coma derived here is in conclusion 
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Figure 3. Panel a) The log(D n )-log(cr) relation for the Virgo, 
Fornax, and Coma samples using the data from LGCT for the 
latter cluster. Panel b) same as panel a with the data of JFK. 
The solid, dotted, and dashed lines give the fits to the Virgo, 
Fornax, and Coma data-points respectively. 

slightly smaller than the results of other DIs. The responsi- 
bility of this difference is likely in the adopted values of the 
central velocity dispersions. Within the uncertainty, how- 
ever, the adopted value is in good agreement with the B- 
band TF relation (Fukugita et al. 1993), with the D n -a re- 
lation of Dressier et al. (1987), with the SN-Ia (Capaccioli 
et al. 1990), and with Bower et al. (1992), who used a com- 
bination of three methods: CM diagram, colour-diameter re- 
lation, and FJ relation. 

Another possible source of error, discussed by LGCT 
and Lucey et al. (1991a), is the surface brightness bias which 
may affect the D n -a relation of the Coma galaxies. Such a 
bias is evident in the LGCT data looking at the residuals of 
the D n -a relation for the galaxies in the central region of 
the Coma cluster with respect to those in the outer parts. 
We checked the influence of this bias on A/icv by selecting 
objects in the core of the cluster and outside, but we did not 
find any significant difference between the two samples. 

3.4 The aperture effects on the central velocity 
dispersion 

Since aperture effects play a major role in determining the 
measured central velocity dispersion, we perform another 
test to better investigate this problem. 

In Fig^ (bottom panel) we plot the difference between 
the velocity dispersion values for the same galaxy measured 
through two apertures of different radius r s < n (small and 
large slit apertures) versus the aperture size increment in 
units of the galaxy effective radius. The data are for the 
Virgo galaxies: the velocity dispersions have been measured 
from slit widths of 2", 4" and 16" and are taken from the 7- 
Sam and Dressier (1984) samples. The solid line represents 
a fit obtained for the galaxies of our Fornax sample. For 



Figure 4. Bottom panel) The open square represent the differ- 
ence between the velocity dispersions of the Virgo galaxies ob- 
tained through the slits of 2 and 4 arcsec respectively. The filled 
square that from 2 and 16 arcsec. The solid line is a fit to the 
galaxies of our Fornax sample (see text). Upper panel) The dif- 
ferences of a from small and large apertures from the catalog of 
Davoust et al. (1985). 



each galaxy we calculate the difference between the central 
velocity dispersion and the luminosity weighted velocity dis- 
persions within increasing fractions of the effective radius. 
The average values for 13 galaxies can be fitted by the rela- 
tion: (To — o ap = 18.3 x log(r ap /r c ) +38.3 with a correlation 
coefficient of 0.97 and a rms of 2.57. 

The upper panel reports the data presented in the cat- 
alog of Davoust et al. (1985). Here the authors listed the 
width of the slit normalized to the corrected face-on diame- 
ter of the galaxies D(0). 

The figure shows that when the difference among the 
apertures is larger than ~ 0.1 r c , the value of a obtained 
through the small aperture is systematically larger than that 
given by the large slit. It should be emphasized that for the 
smaller galaxies this effect is much more evident (the slit 
covers a larger fraction of the effective radius). 

The scatter around zero at smaller A£ = (r; — r„)/r e 
is possibly explained by the 10-15% error of the individ- 
ual measurements, but it may have also another explana- 
tion. According to Ciotti et al. (1996) and Ciotti & Lanzoni 
(1996), the projected central velocity dispersion depends on 
the fraction of effective radius of the galaxy covered by the 
slit width and on the surface brightness distribution of the 
galaxy as whole. The ratio of the central to the virial a as a 
function of r/r e has approximately a gaussian shape, with 
a peak centered between 0.1 0.3r c . The shape and posi- 
tion of the peak depend on the exponent m of the r 1//? ™ law, 
which best fits the light distribution of early-type galaxies 
(Caon et al. 1993). When m is greater than 4, the curve be- 
comes progressively monotonic decreasing. This particular 
behaviour can be responsible of the observed scatter around 
zero in Fig^. In fact, according to the fraction of r c covered 
by the slit one can get different a. Only for the galaxies with 
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m > 4 the difference between the central velocity dispersion 
obtained from the small (a s ) and large (<tj) apertures is al- 
ways greater than zero. Interestingly enough we verified that 
for the Virgo galaxies with m > 4 we always have a s —ai > 0. 
Anyway, taking into account the uncertainties on the mea- 
sured <r, the prediction of this model cannot be presently 
confirmed. 

By applying the aperture correction to the data of 
Dressier (1984) for the Virgo and Coma clusters the result- 
ing D n -<j relation is log(D n ) = 1.531og(<j) — 1.79 and the 
differential distance modulus A/icv= 3.685. The same cor- 
rection applied to the data of the 7-Sam (for Virgo) and to 
our Fornax sample gives log(D n ) = 1.531og(<r) — 1.81 and 
A(j,fv= 0.35. Notice the change of the slope from 1.26 to 
1.53, and the better agreement of the differential distance 
moduli with the average literature values. 

3.5 The error analysis 

As already stated, the limits of the above results in the de- 
termination of the differential distance moduli are in: 1) the 
a priori assumption that the slope of the D n -cr relation is 
the same for the three clusters, 2) the choice of the "rep- 
resentative" sample of objects for the whole cluster, 3) the 
accuracy of the data (in particular ao), and 4) the defini- 
tion of the variables which enter the relation: e.g. ao or the 
average a within a prefixed radius. 

Before turning our attention to the errors associated to 
the use of the D„-a relation, we have checked the assumed 
invariance of the slope dlog(D„)/dlog<7o. To this end we 
have shifted the Fornax and Coma galaxies to the distance 
to the Virgo cluster using the above calculated Alog(D„), 
and fitted the whole distribution again. The result: 

log(A,) = (1-24 ± 0.06) log a - (1.08 ± 0.15), (3) 

with r.m.s. = 0.08, shown in Fig. ^ is in excellent agreement 
with eq. [j]. 

Let us now analize the fitting procedure. Isobe et al. 
(1990) and Feigelson & Babu (1992) developed five meth- 
ods for applying linear regression fits to a bivariate data 
distribution with unknown or negligible measurement errors. 
The advantage of their method is that the uncertainties in 
the slope and intercept of the fitted distribution, and in the 
process of comparative calibration, are evaluated through re- 
sampling procedures: jackknife and bootstrap error analysis. 
The bootstrap error analysis is based "on the distribution 
of slopes and intercepts of a large number of data sets con- 
structed by random sampling of the observed data set with 
replacement". The jackknife error analysis is performed in a 
similar way but with only N synthetic data sets, each con- 
taining (N — 1) points from the original data set, leaving out 
one observation in sequence. 

Following their approach we used the standard linear 
regression analysis (program SLOPES, kindly provided by 
the authors), with the distance dependent variable (log(D„)) 
as the Y variable. This choice has the effect of minimizing 
the residuals in D n . The slopes for the Virgo, Fornax, and 
Coma samples range from 1.16 to 1.37, with different errors 
for the different samples. For the Fornax cluster the errors in 
the intercept (0.69) and in the slope (0.31) are much larger 
than for the other two clusters, due to the large scatter of 
the data and to the small number of galaxies. 




2 2.2 2.4 2.6 

log(tj) [km/s] 



Figure 5. The fit to the full data-set of Virgo, Fornax, and Coma 
galaxies (C 2 D, LGCT, and JFK samples, for a total of 128 ob- 
jects), after shifting Fornax and Coma to the distance to Virgo. 
Symbols are the same as in the previous figures. 

Which is the error that affects the intercept offset ? As- 
suming that the distance to the calibrating cluster is known 
with high accuracy (i.e. that intercept and slope of the fit are 
error free), the error in the relative distance comes approxi- 
mately from the uncertainties in log(D n ) (Ed = Alog(D n )) 
and in log a (E a = A log(cr)) for the new sample of galaxies. 
Since we are looking for the fit that minimizes the residuals 
of the new sample, and we assume that the slope is the same 
as that of the calibrating sample, we can set the maximum 
shift in the intercept equal to ~ * / E D + E% . Being the er- 
rors in both coordinates equal to 0.02 (minimum error), the 
minimum shift of the intercept is ~ 0.03, corresponding to 
an error in A/z of 0.14 mag or 6% in distance. 

A more correct estimate of such errors comes from the 
program CALIB developed by Isobe et al. (1990) and Feigel- 
son & Babu (1992), which predicts the value of the intercept 
(and its uncertainty) in the fit of a new sample of galaxies, 
from the linear regression of a calibrating sample, without 
the assumption that the two slopes are equal. 

In summary this procedure gives A log(D^ v ) = 0.093 ± 
0.015 and Alog(D^ v ) = 0.095 ± 0.06, assuming Virgo 
and Fornax respectively as calibrators, and Alog(D^ v ) = 
0.693 ± 0.013 and Alog(£>£ v ) = 0.690 ± 0.016 for Coma 
and Virgo, A 0.015 uncertainty in log(D n ) corresponds to 
~ 3% error in the relative distance to the clusters. This is 
the minimum uncertainty attainable with such data for the 
relative distance between the two clusters. 

4 THE FUNDAMENTAL PLANE RELATION 

For the Virgo galaxies of the C 2 D sample, the FP relation 
shown in Fig. n writes: 

log(r e ) = (1.26 ± 0.09) log(ff) + 

(0.28 ± 0.01) (^} 8 - (7.31 ±0.28). (4) 
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Figure 6. The FP for the Virgo (open circles), Fornax (filled 
circles), and Coma (open triangles) galaxies. The solid line is our 
eq. ^, the dotted and dashed lines represent the fit of minimum 
scatter for Fornax and Coma respectively, assuming the same 
slope of Virgo. 



The values of the coefficients are derived here through a 
multivariate statistical analysis. The r.m.s. error of the fit is 
0.06, which translates into a 15% uncertainty in the distance 
to a single object, the same obtained with the D n -a relation. 
For comparison JFK found for the Coma cluster: log(r c ) = 
(1.20 ± 0.09) log(a) + (0.35 ± 0.01) (n) 6 - (9.31 ±0.36), with 
a r.m.s. of 0.046. 

Assuming that the slope provided by the Virgo cluster 
is also valid for Fornax and Coma, and using the JFK data 
we found respectively Alog(r^ v ) = 0.10 and Alog(rJ? v ) = 
0.69. Using the seeing corrected parameters of J0rgensen et 
al. (1995a) we found a A/icv = 3.5 mag in good agreement 
with the above uncorrected quantities. 

This means that the FP method gives exactly the same 
relative distance to Virgo obtained through the D n -a rela- 
tion. The same result is obtained for Coma if we adopt the 
FP relation derived by JFK. 

Now, by setting X = 1.26 log(a) + 0.28 (fi) e - 7.31 
and Y — log(r c ), and running again the program SLOPES 
and CALIB, we obtain: Alog(rJf v ) = 0.06 ± 0.01 and 
A log(r^ v ) = 0.07 ± 0.05, using Virgo and Fornax respec- 
tively as the calibrating samples, Alog(rJ? v ) = 0.69 ± 0.01 
and Alog(r? v ) = 0.70T0.02, using instead Virgo and Coma. 
Again the uncertainty in log(r c ) corresponds to a ~ 3% er- 
ror in the relative distance to the clusters. The results of 
these programs are summarized in Table 4. 

If we shift Fornax and Coma to the distance of Virgo 
we obtain for the whole sample of 65 galaxies: log(r c ) = 
1.02-X + 0.07 with r.m.s. = 0.09, comparable to the intrinsic 
scatter of the relation, and corresponding to an uncertainty 
of ~23% in distance. This is a value somewhat larger than 
previously reported in the literature. 

The distance moduli to Fornax and Coma relative to 
Virgo, based on the FP relation, now span the interval 



0.3 < A^ipv < 0.45 mag, and 3.45 < A^cv < 3.50 mag, 
depending on the fitting method that is used. 



5 TESTING THE D n -a AND FP RELATIONS 

The second part of this paper is dedicated to a critical anal- 
ysis of the D n -a and FP relations, taking advantage of our 
well defined (homogeneous and volume-limited) sample of 
galaxies in Virgo and Fornax. 

5.1 The surface brightness bias in the D n —a 

relation for the Virgo and Fornax clusters 

The existence of a surface brightness bias in the D n -a re- 
lation has been claimed by many authors (e.g. Lynden-Bell 
et al. 1988, LGCT, JFK, van Albada et al. 1993). We have 
looked for such a bias using the D n -a relation modified by 
adding a term in surface brightness, as proposed by van Al- 
bada et al. (1993): 

log(A0 = a log (cr) + 6 A/i + cA/i 2 , (5) 

where Afi — {/i) e — 20.75. For our Virgo sample we obtain 
a — 1.356, b = —0.008, and c = 0.002, with a r.m.s. scatter of 
0.06. The multivariate statistical analysis performed on the 
three variables log(D„), log(<r), and A^i, shows that 95% of 
the total variance is enclosed in the correlation between the 
first two variables. For comparison the coefficients found by 
van Albada et al. (1993) with the data of Faber et al. (1989) 
are a = 1.35, 6 = 0.047, and c = -0.019. With our data the 
contribution of the terms in surface brightness is very low. 

We show in Figj^ the distribution of the galaxies of the 
C 2 D and JFK samples in the D n -a-A[i diagram. Again the 
offsets among the clusters distributions are of the same order 
as in the classical D n ~a and FP relations. 

The trend of the residuals Alog(D n ) in the D n -a rela- 
tion with respect to (/x) e is shown in Fig.^ for the galaxies 
of the C 2 D sample. Here we have explored different defini- 
tions of D n by varying the average surface brightness {n) n 
from 19.75 to 22.75 mag arcsec -2 . It is apparent that there 
is no clear trend of the residuals with (/i) c , in agreement 
with LGCT, but not with JFK. The r.m.s. scatter of the 
residuals is similar in all cases and equal to 0.09. 

Following Dressier et al. (1987), who showed that 
D„/2r c cx (I) c ' 8 for r 1/4 galaxies, D'Onofrio et al. (1994; 
their eq. 8 and Fig. 1) derived the relation between D n /2r c 
and (jit)e for objects with r 1//m light profiles. Here m is the 
exponent that provides the best fit to the major axis light 
profiles for the early-type galaxies of the C 2 D sample, and 
ranges from 1 to > 10, increasing systematically with the 
galaxy luminosity. 

Under simple assumptions, the inverse relation {fi) c — 
versus D n /2r e can be written: 

oo 

log(D„/2r c ) = a k (m)[(n) c - (fi)„] k (6) 

k = 

where {fi) n is the average surface brightness within a pre- 
fixed isophote (in general {n)n= 20.75), and the coefficients 
afe(m) are functions of m. 

Writing the observed D n -a relation as: 

log(D„) =a logo- + 6' (7) 
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Figure 7. The D n -cr—Afi relation for the Virgo, Fornax, and 
Coma samples of galaxies. Symbols are the same as in the pre- 
vious figures. The solid line is our eq. ^, the dotted and dashed 
lines represent the fit of minimum scatter for Fornax and Coma, 
assuming the same slope of Virgo. 

and the FP equation as: 

lOg(-Re) = Ox log (J -j- b{fi) c + c, (8) 

the difference between eqsj^ and can be expressed by the 
formula: 

Alog(D„) = [a - a'} loga + [b + ai(m)](^) c + 

OO 

a k (m)[(,u) c - {jj)n] k + const, (9) 

k = 2 

where const=ao(m) + c — b' — 20.75 ai(m). The values of 
the first three coefficients at for m — 4 are approximately: 
a = 1.2 x 1(T 3 , ai = -2.7 x 10~\ a 3 = -1.5 x 10" 2 , 
vanishing for large values of k. 

Looking at the values of the slope and of the intercept of 
the D n -a and FP relations, one sees that all the coefficients 
in eq. ^ are very small (e.g. [a — a'] ~ 0.05, assuming respec- 
tively the slopes in eq. ^ and eq. [|, [b + Oi(4)] = 0.007, and 
so on) for each value of m. For this reason the dependence 
on (/i) c of the residuals of the D n -a is almost completely 
cancelled out, and this is confirmed also in our Fig. H. 

5.2 The behaviour of the residuals in the D n —a 
and the Fundamental Plane relations 

A test of the behaviour of the residuals of the D n -a and 
FP relations is presented in Fig. || where the residuals 
Alog(D„) and Alog(r e ) for the C D galaxies are plotted 
versus the galaxy luminosity Mb, the ellipticity at the effec- 
tive radius e e , the normalized coefficient of the Fourier 
expansion of the residual of the isophotes with respect to 
the best fitting ellipse (Bender & Mollenhoff 1987), the log- 
arithm of the exponent m a of the r 1,/m law that best fits the 
major axis light profiles of the galaxies, the logarithm of the 



maximum observed rotation velocity along the major axis 
of the galaxies Vm, and the logarithm of the sum V m + ao- 
No clear trend of the residuals is seen in almost all cases, 
with the possible exception of the galaxies with large values 
of m a and low maximum rotation velocities Vm- galaxies 
with Vm < 60 km s _1 have preferably negative residuals, 
the opposite being true for the faster rotators. Unfortunately 
this result rests on a small number of objects. The Vm values 
for the Fornax galaxies have been derived from the rotation 
curves measured by D'Onofrio et al. (1995), while for Virgo 
we used the sub-sample of 12 objects with the best observed 
rotation curves and velocity dispersion profiles available in 
the literature. 

Notably, the small asymmetry in the residuals disap- 
pears if we consider the quantity V m + a 2 , which is more 
closely related than V m to the total kinetic energy. 

The dependence of the residuals on the kinematical 
properties of the galaxies is more clearly seen in Fig. |Io[ 
where Alog(D n ) is plotted versus the ratio V/a between 
the ordered and random motions, and versus the anisotropy 
parameter (V/a)* (cf. Davies et al. 1983), which is a measure 
of the departure of an object from the rotational support. 
In the lower plots V/a is computed using the highest values 
assumed by the rotational velocity and the velocity disper- 
sion in the observed range, while in the upper plots it is the 
ratio of the values of V and a measured at the galactocentric 
distance 0.3a e , where a c is the effective semi-major axis of 
the galaxies. 

There is possibly a trend in the residuals: we find 
Alog(Ai) = 0.091og(V m /<7) + 0.06 with a r.m.s. equal to 
0.1 and a correlation coefficient of 0.34. ^From the t-Student 
distribution the (null) hypothesis of no correlation between 
the two variables can be rejected at the 95% confidence 
level. This mild correlation, however, is not completely at- 
tributable to the presence of the SO galaxies, but rather is 
due to the different distribution of the galaxies dominated 
by anisotropic pressure (log(Vm/o")* < —0.25) with respect 
to those supported by rotation (\og(Vm/a)* > —0.25). In 
fact, if we eliminate the SO galaxies, a slight trend is still 
apparent: Alog(D n ) = 0.06 log(Vm/o") + 0.02, with a corre- 
lation coefficient equals to 0.21. The trend in the residuals 
disappears in the correlation with log(V/o")o.3a e when the 
SO are eliminated. At large radii only the SO galaxies, that 
are largely supported by rotation, have positive residuals. In 
the upper plot only few Virgo galaxies have rotation curves 
and velocity dispersion profiles extended to 0.3a c . 

In the above discussion, we did not take into account 
projection effects in determining V m , since this would require 
the knowledge of the intrinsic structure of each galaxy. This 
is of course incorrect. However, Binney (1978, 1980) showed 
that for oblate ellipticals with Hubble luminosity profiles, 
the ratio {V 2 )/{a 2 } can be approximated to better than 20% 
by V 2 q /a 2 (where V 2 q is the true maximum equatorial rota- 
tion velocity and a c is the velocity dispersion measured at 
the centre of the galaxy), and that V 2 q /a 2 — (16/iv 2 )Vo /ao, 
(where Vq and ao are the projected values for an edge-on 
oblate galaxy) . Since we can assume a random orientation of 
the rotation axes with respect to the line of sight, the mean 
projection of V 2 along this line is V m = 2/3V 2 , and we can 
believe that, to a first order approximation, the use of Vm 
and ao should not produce a systematic error of (V 2 )/(a 2 ). 

Projection effects along the line of sight are also impor- 
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Figure 8. Residuals of the D n -a relation versus the average surface brightness (fi) e for the galaxies of the C 2 D sample. Different choices 
of (yu)n are used in each box (see text). Open and filled symbols represent Virgo and Fornax objects respectively. Circles and triangles 
are E and SO galaxies respectively. The dashed lines are best fits to the plotted distributions. 



tant. In this case rotation velocities are generally underesti- 
mated (by a factor 1.5-2), while velocity dispersions remain 
almost unchanged (Busarello et al. 1992). We have not tried 
here any deprojection since this would require the knowl- 
edge of the rotation axis of the galaxy with respect to the 
line of sight. The measured V/a ratios are therefore lower 
limits. 

Unfortunately the small number of data does not permit 
a full analysis of the effects of rotation on the D n -a and FP 
relations. We can only state that our result is along the 
same line of that of Saglia et al. (1993), who claim that 
the presence of disk components inside ellipticals produce a 
scatter in the FP, and is in agreement with that of Prugniel 
& Simien (1994), who found the same correlation of V/a 
with the residuals of the FP, using a sample of galaxies which 
is not volume-limited. 

Indeed, despite the poor correlation, it seems that the 
kinematical structure of the galaxies plays an important role 
in the D n -a and the FP relations, in particular for what con- 
cerns the thickness. The amount of the detected deviations 
is in agreement with the prediction of the Virial theorem, 
when the contribution of rotation to the global equilibrium 
of elliptical galaxies is taken into account. 



Fig. |llj shows that the rotationally supported galaxies 
provide a mean relative distance modulus for the Fornax 
cluster with respect to Virgo which is lower (A/ipv = 0.2 
mag) than that obtained from the galaxies dominated by 
velocity anisotropy (A/xfv = 1.0 mag). Actually, this result 
is uncertain due to the large scatter of the Fornax galaxies 
and the small number of objects, and should be considered 
here only as an example of the bias that can be introduced by 
the rotation effects. If confirmed by future data, this source 
of bias should be taken into account in the D n -a and the 
FP relations when used as DIs of high redshift clusters. The 
misclassification of SO galaxies would in fact introduce a 
systematic bias. 

Can we correct the D n -a for such effect ? A simple way 
could be the following. Let us write that D n cx {a 2 ) a , which 
means that the diameter of a galaxy is strictly related to the 
central velocity dispersion. If the contribution of the rotation 
cannot be neglected, we can write that D n ' <x (a 2 + f3V 2 ) a , 
where j3 parametrizes the degree of rotational support. Tak- 
ing the logarithm of these quantities and assuming that 
a' — a ± Aa, one obtains: 

Alog(.D n ) = (q± Aa)log(l + /3fc 2 ), (10) 
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Figure 9. The residuals Alog(D n ) and Alog(r e ) of the D n -cr and FP relations with respect to the absolute magnitude Mb, the 
cllipticity at r c , log(ra a ), the boxy/disky parameter (0,4/0) * 100, log(V^ + cr 2 ), and the maximum observed rotation velocity along the 
major axis of the galaxies log(V m ). The Virgo and Fornax galaxies are indicated by open and filled symbols respectively. 
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where k — V/a. 

Unfortunately we do not know the correct values for the 
parameters a and 8. The rotational support, parametrized 
by 8, is different for each galaxy, but a mean values for early- 
type galaxies could be within 0.5 and 1. The same interval 
is probably spanned by a (for the classical D n -a relation 
a ~ 0.6). Using the approximation a — 1 and B — 0.5, one 
obtains the correction curves plotted in Fig. |l2|. 

After applying this correction, the trend of the resid- 
uals of the D n -a relation disappears, but the slope of the 
fitted distribution increases from 1.26 to 1.37, thus affecting 
the distance determinations of the clusters. Unfortunately 
the r.m.s. of the fitted distribution does not decrease signif- 
icantly. 

By considering only the sample of galaxies possessing 
both oo and V m , we find that these two variables are not 
mutally correlated and that log(Vm/<7)* does not correlate 
with Mb, while it correlates with (/n) e (as in Wyse & Jones 
1984). It is also interesting to note that the residuals of 
the D n -a and FP relations, Alog(D n ) and Alog(r c ), are 
mutually correlated. We find Alog(r e ) = 0.94A log(D„) + 
0.05 with an r.m.s. of 0.036. This is an obvious consequence 



of the Dressier et al. relation D n /2r c oc (/) c °' • The two 
relations are therefore almost equivalent when used as DIs 
for nearby clusters. 



5.3 Further test on the D n —a relation 

Do variations in the definition of D n (i.e. in the surface 
brightness level (fi) n ) affect the D n -a? We performed an 
experiment with the Virgo and Fornax C 2 D sample by vary- 
ing the value of D n , i.e. by assuming that {fj) n changes from 
= 18.75 to 22.75 mag in steps of 1 mag. In the case of Virgo 
the slope is found to have small variations, from 1.21 to 1.31, 
while for Fornax the variation is larger, and range from 1.19 
to 1.37. This is due to the larger spread in the central ve- 
locity dispersions. The r.m.s. of the fit varies from 0.07 to 
0.11 for Virgo and from 0.13 to 0.16 for Fornax. The best 
fits occur at hb = 19.75 and 20.75 mag arcsec -2 . All these 
expressions for the D n -a relation provide approximately the 
same A log(D n ) (it varies from 0.09 to 0.11), that is the same 
Afi. 

In another experiment we have acted on the definition 
of a, adopting once the value of the velocity dispersion mea- 
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Figure 11. The D n -o relation for the Virgo and Fornax clus- 
ters. The galaxies of Fornax are plotted with different symbols 
according to the value of log(V7 cr) max . Two Fornax galaxies with 
<jo < 100 km s _1 have been added here to increase the number 
datapoints. The larger scatter of the Fornax data with respect to 
Virgo is more evident in this plot. 
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Figure 12. The solid line is the correction curve for the D n 
diameters according to eq. |Io[ with a = 1 and j3 = 0.5. The 
dashed lines give the curves for Act = ±0.25 and the same /3. 



sured at the centre, ao, and another time the average within 
a given radial range. The motivation for this experiment 
comes from the Virial Theorem, since we want to look for 
the value of a that best represents (a). 

We decided to average a within circles of radius r n = 
D n /2 and 0.3a c , where a c is the effective semi-major axis. 
In each case we took plain averages, weighted means for the 



surface brightness, and means weighted for the rotational 
contribution. Unfortunately, this test was possible only for 
the 15 galaxies of the Fornax cluster which have accurate 
extended rotation curves, velocity dispersion profiles, and 
luminosity profiles. We found that the slopes of the D n -a 
relation obtained in each case vary in the interval 1.34-1.62 
(to be compared with the slope of 1.26 achieved using the 
central velocity dispersion) , that is it changes by an amount 
larger that the error of the fit, and increases with the radius 
out to which the average is computed. Approximately equal 
r.m.s. fits are obtained using the central velocity dispersions 
ao and the value averaged within 0.3a c . Slightly larger errors 
are derived with the means within r n . 

These results are an indication that early-type galaxies 
deviate from the pure homology and that the D n -a and the 
FP relations have a slope and an intrinsic scatter connected 
to the spatial and dynamical structure of the galaxies. 

One might try to use the quantity (V^ + a 2 ) instead of a 
in order to account for the above mentioned bias introduced 
by the different kinematical properties of the galaxies. The 
result of such an experiment is that the slope of the relation 
changes from ~ 1.23 to ~ 1.37, but the r.m.s. of the fit does 
not show a significant improvement. 



6 THE log(m)-log(r c ) RELATION 

We finally present another method to derive the relative dis- 
tance to the Virgo and Fornax clusters. This method is based 
on the log(m)— log(r e ) relation (Caon et al. 1993, Young & 
Currie 1994) and is shown in FigJ^j. Here m a is the ex- 
ponent of the r 1 /" 1 Sersic' law best fitting the major axis 
surface brightness distribution of the galaxies, and R a is the 
effective radius in kpc. Although the scatter in this relation 
is quite large (the r.m.s. of the fit is = 0.28 and the corre- 
lation coefficient is 0.7), the intercept offset Alogi? c turns 
out to be 0.089, corresponding to A/xfv = 0.44 mag, in 
excellent agreement with the D n -a and FP relations. Such 
agreement is remarkable although probably fortuitus given 
the large scatter of the data in this relation. 



7 CONCLUSION 

The numerical results obtained in this work fitting the D„-a 
and FP relations are summarized in Table 4. /.From these 
data we derive Afipv = 0.45 mag and A(j,cv ~ 3.45 mag 
respectively for the relative distance moduli to Fornax and 
Coma with respect to the Virgo cluster. Taking into account 
the corrections for the aperture effects A/^py is closer to 
0.35 mag and A/^cv to 3.55 mag. Our distance to Fornax 
is slightly larger than derived with the other DIs, but it 
is nearly equal to that derived with the BCM method (see 
Table 1). 

The inner ring diameters of spirals, the globular cluster 
luminosity function (GCLF), the surface brightness fluctu- 
ations (SBF), the L—a-fj, relation, the IR Tully-Fisher re- 
lation (IR-TF), and the scale length of dE galaxies, are in 
open contrast with this determination, while the planetary 
nebulae luminosity function (PNLF), and the SNe-Ia are 
consistent. 

For the Coma cluster we derived a distance smaller than 
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Table 4 



Method 


a 


b 


c 


rms 


A/iFV 


A^cv 


Notes 


Dn-cr 


L26 ±0.06 


-1.11 ±0.14 




0.06 


0.55 




1 


D n -a 


1.31 ±0.07 


-1.24 ±0.16 




0.06 


0.45 


3.40 


2 




1.24 ±0.06 


-1.08 ±0.15 




0.08 


0.45 


3.45 


3 


Dn-a 


1.26 ±0.06 


-1.11 ±0.14 




0.06 


0.46 


3.46 


1 




1.37 ± 0.31 


— 1.46 ± 0.69 




0.13 


0.47 




5 


Dn-a 


1.16 ±0.08 


-1.59 ±0.19 




0.08 




3.45 


6 


L—a—[i 


1.35 ±0.08 


-0.008 ±0.178 


0.002 ± 0.016 


0.08 


0.45 


3.45 


7 


FP 


1.26 ±0.09 


0.28 ±0.01 


-7.31 ±0.28 


0.06 


0.45 


3.45 


8 


PP 


1.01 ±0.03 


0.07 ±0.04 




0.06 


0.30 


3.45 


9 


FP 


0.85 ±0.10 


0.23 ±0.13 




0.13 


0.35 




10 


FP 


1.15 ±0.04 


-0.84 ±0.06 




0.06 




3.50 


11 


log(m)-log(r e ) 


1.16 ±0.25 


-0.160 ±0.14 




0.28 


0.44 




12 



Notes: 1) Linear least square fit for the Virgo galaxies and fit of minimum scatter for Fornax adopting the same slope of Virgo. 2) Fit 
to the Virgo galaxies taking into account the errors on both coordinates. 3) Fit of minimum scatter for Fornax and Coma using the fit 
of the combined sample of 128 galaxies. 4) Output of program CALIB using Virgo as calibrator. 5) Output of program CALIB using 
Fornax as calibrator. 6) Output of program CALIB using Coma as calibrator. 7) Coefficients of the L-a-fi relation from the 
multivariate statistical analysis. Fit of minimum scatter for Fornax and Coma, adopting the same slope. 8) Coefficients of the FP 
relation from the multivariate statistical analysis. Fit of minimum scatter for Fornax and Coma, adopting the same slope. 9) Output of 
program CALIB using Virgo as calibrator. 10) Output of program CALIB using Fornax as calibrator. 11) Output of program CALIB 
using Coma as calibrator. 12) Coefficients of the fit for the Virgo galaxies. Fit of minimum scatter for Fornax, adopting the same slope. 
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Figure 13. The log(m)-log(r e ) relation for the Virgo and Fornax 
galaxies. R c is the effective radius of the galaxies in kpc. The solid 
line gives the best fit for the Virgo galaxies. The dashed line is 
the fit of minimum scatter for the Fornax sample. 



given by the other DIs, but now a better agreement exists 
with the TF relation. 

Probably our estimates can be improved when new ac- 
curate and homogeneous central velocity dispersions and ro- 
tation velocities will be available for a complete sample of 
Virgo and Fornax early-type galaxies. Actually, the large 
intrinsic scatter we found for the D n -cr and FP relation 
of Fornax is at variance with previous determinations. We 
believe that the problem is in the measurements of a: the 
large spread of our central velocity dispersions being caused 



by projection effects, a anomalies, and the presence of inner 
sub-components such as disks or central black holes. In this 
context, the large uncertainty on Afi resides probably in the 
different procedures used for measuring the velocity disper- 
sions and in the corrections adopted for taking into account 
the effects of the aperture sizes. 

By applying an average aperture correction to the cen- 
tral velocity dispersions, we obtain differential distance mod- 
uli in better agreement with the literature data, and we find 
that the tilt of the D n -cr relation increases. 

Looking at the data presented in this paper we conclude 
that the D n -o and the FP relation give substantially the 
same result. An advantage of the D n -a relation over the FP 
is that D n can be measured with high photometric accuracy, 
while the effective radius r c and {fi) e may depend on the 
actual shape of the luminosity profiles and on the error in 
the total luminosity of the galaxies. The errors in r c and 
{/x) e are in fact usually correlated (Capaccioli et al. 1992). 

At variance with van Albada et al. (1993), we do not 
find any need for adding a surface brightness correction to 
the D n -u relation, at least for the galaxies of the C 2 D sam- 
ple. The inverse (D n /A e )-{fj,) c relation of Dressier et al. 
(1987), coupled with the observed D n -a and FP relations, 
has indicated that the residuals of the D„-ir relation depend 
only slightly on {/i) c , since all the coefficients in eq. |^ are 
approximately null. 

The residuals in both the D n -a and FP relations are 
found to be independent of the intrinsic luminosity of the 
galaxies, the ellipticity, the mean surface brightness and the 
boxiness/diskiness parameter 04. A trend of the residuals, 
on the other hand, is present with the exponent m of the 
r i/m g ers j c ' s j aW] w jth the maximum rotation velocity of the 
galaxies and with the Kn/cr ratio. The last result is in agree- 
ment with Busarello et al. (1992) and Prugniel & Simien 
(1994), who first pointed out the contribution of rotation to 
the global equilibrium of elliptical galaxies. 

The correlation of the residuals of the D n —cr relation 
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with the (V/a) ratio, with the anisotropy parameter (V/a)*, 
and with the exponent m a of the major axis light profiles of 
the galaxies are all indications of the non-homologous nature 
of the early- type galaxies. 

Further improvements call for more accurate kinemat- 
ical measurements, for a critical estimate of the different 
methods of data reductions, for a standard correction of the 
aperture effects, and for a deeper analysis of the kinematics 
of each individual galaxy. Even at this stage, we are probably 
legitimated to say that most of the observed deviations from 
the FP are likely due to departures from the pure homology 
in the density distribution and in the dynamical structure 
of the early- type galaxies. 

An important step for the use of the D n ~a and the FP 
relations as DIs will be that of defining a "representative 
sample" of galaxies, of using standard selection criteria, and 
data with the same degree of homogeneity and accuracy, as 
well as variables properly and univocally determined. 

In closing we want to remark that, while our analysis 
disagrees with Saglia's et al. (1993) for what concerns the 
dependence of the residuals on the ellipticity and on a,4, 
nonetheless we back their conclusion: the kinematical struc- 
ture of the early-type galaxies, if not properly taken into 
account, may produce biased cluster distances, when the 
D n -a and FP relations are used as DIs without care. 
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